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Abstract 

We describe the morphology of the universal function for the critical 
(cubic) circle map at the golden mean, paying particular attention to the 
birth of inflection points and their reproduction. In this way one can fully 
understand its intricacies. 

Introduction 

At the critical point, the circle map has a cubic order of inflection. Near the 
golden mean winding number the self-reproducing (renormalization group) prop- 
erties of the Fibonacci mode-locked steps [1] leads to the notion of a universal 
circle map function g{x), which has received a great deal of study in the litera- 
ture [2, 3, 4]. It obeys the compatible [5] pair of functional equations, 

g(g(x)) = g(a 2 x)/a, (1) 

g(g(x)/a 2 ) = g(ax)/a 2 . (2) 

The universal circle map constant a = —1.288575.. [2] can be determined to 
extremely high accuracy [6] by examining the equations over a restricted range 
of x. In this paper however we focus on the behaviour of g(x) over large x, in 
order to make sense of its morphology. Some years ago we tried a similar analysis 
for the period-doubling universal function [7]; our aim here is to comprehend 
how the ever-increasing number of inflection points are spawned as we go out 
to ever larger We shall see that they can be discovered directly by a close 
look at the functional equations and their derivatives, starting from previously 
known series of inflection points. 

In the next section we quickly recapitulate the origin and character of the 
functional equations and derive the principal features, including the first family 
of inflection points. In the last section we locate the higher families of inflection 
points of g and explain how successive ones are derived theoretically. In this way 
one can better appreciate the morphology of the universal circle map function 
and comprehend its complicated structure. 

Origin and Main features 

The approach to the irrational golden mean winding number in nonlinear maps 
of an angular variable 9 is normally made [1] via a succession of Fibonnaci 
number ratios: F n —i/F n asn^ oo. In particular for the critical circle map, 

0^0' =/(0,fi) = O + (9 - sin(27r60 /27r, (3) 
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with a cubic inflection point at the origin, it is conventional to work out the 
superstable values of the driving frequencies fi n which lead to the rational 

F n -i/F n , 

[/] F - 1 (0,O„)-F„=0. (4) 

(Above we are using the notation, [f] N (x) to denote the AT- fold composition, 
/(/(■•• f(x) • • •))•) Correspondingly, 

S v ' 

N 

Ae n = [ff—(o,n n )-F n ^ (5) 

is the nearest fixed point to the origin. The universal Feigenbaum constants S, a 
of the map are then obtained [2] as 

6 = lim AQ„/Af2 n _i = -2.833612..; Afi n = fi n -fi„_i, (6) 

n—>oo 

a = lim A6»„/A0„_i = -1.288575.. (7) 

n^oo 

Parenthetically, we note the numerical value fioc = 0.6066610635.. 
By considering the set of functions 

MO) =a n ([/] F "(O n+r ,0/a") -F n _i) , (8) 

in the limit of large n and taking the limit as r —> oo, one may establish that the 
limiting (Shenker) function <fi, evaluated at Ooo, obeys two compatible functional 
equations, 

4,((j)(9)/a) = 4>{a6)/a 2 and (9) 

<f>(auf>(0)) = <t>{a 2 6)/a. (10) 

Numerical approximations to (8) indicate that the function <p(x) (called f(x) 
by Shenker) is a monotonically increasing function of its argument; therefore 
the alternative function g(x) — a<fi(x) decreases monotonically with x and one 
readily establishes that it satisfies the universal equations (1,2). At the same 
time one may easily show that the scaled function G(x) = g(px)/p obeys the 
self-same equation pair (9,10); thus one is at liberty to set the normalization 
scale at will. We shall fix the unique zero of g(x) at x = 1, that is g(l) = 0, 
and thereby find the value g(0) = A = 1.2452.. to be the intercept at the origin. 
Because of the cubic nature of inflection, .g'(0) = g"(0) = 0. 

Plots of g(x) in the ranges -8 to -4, -4 to +4, +4 to +8 are respectively given 
in Figures 1, 2 and 3. These plots are approximate and they were obtained from 
equation (8) by going to high order (n = 14) at the accumulation point and 
scaling appropriately in x, as indicated above. (The range -8 to 8 was chosen 
so that we have sufficiently many inflection points to make the investigation 
of the morphology worthwhile, but not too many to confuse the subsequent 
discussion.) Two significant points may be noted from these plots: 
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• they exhibit a scries of inflection points which are self-similar and which 
proliferate on ever smaller (and larger) scales as we go out in x, 

• the function behaves asymptotically as g(x) ~ ax. 

The last fact is readily verified from the initial functional equations (1) and 
(2); our goal is to understand the first fact, namely the origin and families of 
inflection points. Before doing so let us derive a number of useful facts about g 
at particular locations, some of which are required later. By evaluating (1) at 
x = 1/a 2 and (2) at x = 1/a, one readily discovers that 

g(l/a 2 ) = 1, g(l/a) = a 2 ; 

the former also follows from (2) evaluated at x — 1/a 2 . Taking (1) and (2) at 
x = 1, leads to 

g(a 2 ) = aX, g{ a ) = ce 2 X, 

the former following from (la) evaluated at x — 1/a. Finally, we note that 
working out (1) and (2) at x = 0, gives 

g(X) = X/a, g(X/a 2 ) = A/a 2 . 

The latter condition demonstrates that X/a 2 is the (single) fixed point of the 
universal circle function! Further relations between g(Xa M ) can be found by 
substituting suitable x = a N in cqs. (1), but we will not need them below. 
As the next step, differentiate equations (1,2) to obtain the related pair, 

g'(x)g' (g(x)) = ag'(a 2 x) (11) 

9 ! (x)g' (g(x)/a 2 ) = ag'(ax). (12) 

Remember that where g has zero first derivative, its second derivative also 
vanishes (because the map is cubic), so equations (11,12) provide information 
about the inflection points. A number of interesting values of the derivatives 
can be obtained immediately. Setting x — 1/a 2 in (11) and x = 1/a in (12) 
yields 

g'(l/a 2 ) = a, g'(l/a) = a 

and one can check the correctness of these values from the graphs. Also, taking 
the limit as x —> in (11,12), one deduces that 

ff'(A)=a 5 , g'(X/a 2 )=a\ 

because the inflection point is cubic. By similar means one may establish rela- 
tions between derivatives of g at various points x = a M X which provide useful 
checks on numerical work. 

But one can go further through the following observation: if £ corresponds 
to an inflection point, (11) and (12) ensure that a£ and a 2 £ are also inflec- 
tion points; by induction one generates an entire sequence of such points, viz. 
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a n £, n = 0, 1, 2, However equation (11) at x = 1 and (12) at x = 0, inform 

us that 

sV) = o, </(«) - o, 

so we conclude that = a™, n = 1,2,3, .. . represent a family of inflection 
points, seeded by £i = a; for our purposes we shall regard them as the primary 
or parent series. These points can be picked out in Table 1, where we have listed 
all the inflection points that occur between x = — 8 to +8 and they can also be 
spotted in the figures as marked dots. This table contains many more inflection 
points and the question is: how are they seeded? The answer will be provided 
in the next section. 

Families of inflection points 

Examination of the figures shows that not all inflection points belong to the 
primary sequence a N . The first non-standard inflection point occurs at x ~ 
2.25 and to see where this comes from we examine first the seed values xn = 
g^ 1 (a N ), N = 1,2, ... . These may be read off from the graph of g and are 
tabulated in Table 2. At those locations, 

9' (g(x N ))=g'(a N ) = 0, 

and from (11) we are assured that the set a M XN, M = 2,3, .. . will generate 
a family of inflection points, some of which may not be new. [For example we 
note again that the location x 2 = 1 /a simply generates the primary or parent 
family] The first daughter family, which is new, is given by the secondary series, 
a 2 xi, a 3 xi, . . .; another daughter family, which is also new, corresponds to the 
sequence, ax 3 , a 2 x 3 , . . ., etc. and it starts off with M = 1, by virtue of equation 
(12). However the family a n X4 is not new because one may readily establish 
that X4 — axi by inserting x = x\ in equation (12). New secondary families 
a"a;5, ct n x§, . . . etc. begin with n = 1, again via equation (12). 

While the parent and daughter sequences lead to two generations of inflection 
points, some of which are listed in Table 1, this does not exhaust the class of 
such points. Proceeding in a similar fashion to above, define the granddaughter 
sequence, 

x MN = g- 1 (a M x N ) = g- 1 (a M g- 1 (a N )) . (13) 

Many of these points are not new; by manipulating (8a) and (8b) appropriately, 
one may prove that some daughters are identical with or related to their aunts. 
For instance, 

ax 2 i — x 3 , ax 31 — x 6 , ax 2 5=x r , x 15 — x 6 , X4 3 = ax 6 ,... 
As well, some of these granddaughters are interrelated; for instance 

axi 3 = Xn, ax w = x 5 i, . . . 
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Only the independent a n XMN occurring between -8 and +8 are stated in Table 
1. 

Practically all the remaining missing points of inflection in that range arise 
from the next generation, xlmn = g 1 (cv L xmn) ■ Again we have only listed 
those that spawn new series. The final 'missing' point within -8 to 8, at x ~ 
—7.45, is noteworthy because it belongs to the 'fourth' generation. By following 
through this procedure one may track all the inflection points over any desired 
range, with each family being generated from the previous one, though many are 
identified with previous members; for example the point a 2 X2Ai is the same as 
its aunt ax^i- Within the restricted range -8 to 8, there are no other inflection 
points beyond the ones noted in the table. An interesting observation is that the 
number of inflection points between a n and a n+2 equals the Fibonnaci number 
F n ; while we can readily show that this holds for subcritical maps where a is 
the inverse of the golden mean ratio, an exact proof eludes us for critical cubic 
maps. 

All told we can see how the functional equations reproduce their structure 
in a self-similar way, becoming ever more intricate and extended at one and the 
same time: for instance, the region between -7.7 and -7.4 looks very much like 
the region between -2 and 2. This feature happens everywhere and is not unlike 
the morphology of the period-doubling universal function. We expect that the 
cubic nature of inflection is not vital for the validity of this conclusion and that 
it is common to all scaling functional equations of the type [g] A '(a;) oc g(px). 
The exquisite self-similarity on every scale as we go out in x means that it is 
effectively impossible to make an analytical approximation to g over a wide 
range of x; but of course over a small, limited range it is always possible to 
do so and indeed this is a sensible way to compute the scaling constant a, by 
focussing on the region x = to 1 say. 
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Table 1: Positions of successive inflection points of g{x) to two decimal points, 
between x = — 8 & +8. ID is the interpretation of their origin in the recursive 
notation Xmn... = <7 _1 (a M xn...)- 



X 


9 


ID 


X 


9 


ID 


X 


9 


ID 


-7.99 


10 


60 


a' ' x\ 


-4 


80 


6 


43 


ot b x\ 


3 


73 


-4 


58 


4 

a x\ 


-7.76 


10 


10 


a 2 x 6 i 


-4 


52 


5 


92 


a 3 x\3 


4 


12 


-5 


12 


axe 


-7.65 


9 


91 


a 2 xi6i 


-4 


39 


5 


71 


ax 33 


4 


58 


-5 


46 


a 6 


-7.56 


9 


83 


a 5 x 13 


-4 


27 


5 


62 


a 3 x 3 


5 


12 


-6 


42 


a 2 x 5 


-7.45 


9 


63 


ax 33 i 3 


-3 


98 


5 


26 


ax$ 


5 


27 


-6 


77 


axi4 3 


-7.40 


9 


56 




-3 


55 


4 


99 


a 5 


5 


48 


-6 


91 


4 

a x 3 


-7.35 


9 


52 


a 3 x 33 


-2 


90 


3 


97 


0?X\ 


5 


66 


-7 


22 


a 2 x 33 


-7.26 


9 


39 


a 2 x 233 


-2 


58 


3 


42 


OLX 3 


5 


74 


-7 


37 


a 2 x 133 


-7.09 


9 


30 


a 5 x 3 


-2 


14 


3 


20 


a 3 


5 


83 


-7 


42 


a 4 xi3 


-6.87 


8 


95 


a 2 X43 


-1 


29 


2 


08 


a 


5 


99 


-7 


64 


axQi 


-6.78 


8 


78 


a 2 xi4 3 







1 


25 


1 


6 


19 


-7 


78 


a 6 xi 


-6.62 


8 


69 


a 3 X5 


1 


66 


-1 


62 


a 2 


6 


52 


-8 


29 


a 3 xi6 


-6.34 


8 


32 


OLX-J 


2 


25 


-2 


67 


a 2 x\ 


6 


67 


-8 


53 


ax 36 


-5.90 


8 


06 


a 7 


2 


76 


-3 


08 


a* 


6 


85 


-8 


65 


a 3 XQ 


-5.31 


7 


06 


a 2 XQ 


3 


30 


-4 


07 


a x 3 


7 


20 


-9 


09 


axs 


-5.07 


6 


63 


a 2 x 16 


3 


52 


-4 


41 


a 2 x 13 


7 


60 


-9 


40 


a 8 



Table 2: The primary family of inflection points between -8 & +8 and their 
corresponding inverses (to two decimal points). 



N 


1 


2 


3 


4 


5 


6 


7 


8 


a N 


-1.29 


1.66 


-2.14 


2.76 


-3.55 


4.58 


-5.90 


7.60 


xn = 9~ L (a N ) 


1.35 


-0.78 


1.99 


-1.74 


3.09 


-3.20 


4.92 


-5.58 



Figure Captions 



6 



Figure 1. g(x) from x = —8 to x = —4. This range contains 20 inflection points. 
Figure 2. g(x) from x = — 4 to x = +4. This range contains f 3 inflection points 
(only 8 points between -3 to 3). 

Figure 3. g(x) from x = +4 to x = +8. This range contains 15 inflection points. 



7 




8 



Figure 1: R Delbourgo 




Figure 2: R Delbourgo 
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Figure 3: R Delbourgo 



